Effect of Uniaxial Strain on Ferromagnetic Instability and Formation of Localized 

Magnetic States on Adatoms in Graphene 
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We investigate the effect of an applied uniaxial strain on the ferromagnetic instability due to long- 
range Coulomb interaction between Dirac fermions in graphene. In case of undeformed graphene 
the ferromagnetic exchange instability occurs at sufficiently strong interaction within the Hartree- 
Fock approximation. In this work we show that using the same theoretical framework but with an 
additional applied uniaxial strain, the transition can occur for much weaker interaction, within the 
range in suspended graphene. We also study the consequence of strain on the formation of localized 
magnetic states on adatoms in graphene. We systematically analyze the interplay between the 
anisotropic (strain- induced) nature of the Dirac fermions in graphene, on- site Hubbard interaction 
at the impurity and the hybridization between the graphene and impurity electrons. The polarization 
of the electrons in the localized orbital is numerically calculated within the mean- field self- consistent 
scheme. We obtain complete phase diagram containing non- magnetic as well as magnetic regions 
and our results can find prospective application in the field of carbon- based spintronics. 

PACS numbers: 71.23.An, 73.20.Hb, 75.10.Lp, 75.20.Hr, 75.70.Ak 



I. INTRODUCTION 

Magnetism in organic (carbon- based) materials 1 ^ is 
not only an interesting physical property for studying 
from the fundamental point of view but it is also 
crucial in realizing carbon- based room- temperature 
spintronics-, where one utilizes the charge as well as 
spin degrees of freedom of electrons. Carbon is long 
lasting as well as abundant, and devices made out of 
carbon are expected to be advantageous in many ways 
and inexpensive. Moreover, there are many advantages 
of using organic instead of inorganic substances to make 
spintronics devices. The tunability of electronic prop- 
erties, mechanical flexibility and weak spin scattering 
mechanisms in carbon- based materials could improve 
the reliability of this emergent technology. 

Carbon exists in many allotropic forms. Its two di- 
mensional (2D) form, graphene, consists of carbon atoms 
arranged on a honeycomb lattice and in 2004, it was 
isolated for the first time from graphite in a controlled 
manner 4 . This naturally occurring crystalline pure 
material provides an avenue of many intriguing quantum 
phenomena due to its low dimensional structure. And 
its exceptional mechanical 5 -, electronic^, transports and 
optical- properties has made it a promising material 
for various practical applications. Recently electrical 
spin current injection and detection in graphene was 
established up to room temperature 9 . And owing to 
the absence of nuclear magnetic moments in carbon and 
negligible spin- orbit interaction 10 , graphene is a also 
good candidate for spintronics. 

But the lack of finite electronic band gapii has been 
a hindrance to use graphene for technological purposes. 
Since it has outstanding mechanical properties and is ca- 
pable of sustaining huge atomic distortions 5 , there have 



been proposals to use uniaxial strain as a tool in order 
to open a band- gap and tailor the electronic properties 
of graphene. Pereira et alX^ analyzed such concepts 
within the standard tight- binding (TB) approximation 
and effect of strain was also studied experimentally 
using Raman spectroscopy^ or by using metal cluster 
super-lattices and patterned modifications^ 4 -. Although 
uniaxial strain can not generate the required bulk band- 
gap but there are other important consequences of it on 
graphene besides adjusting the electronic properties 15 . 
For instance, deformations due to strain can alter its 
optical properties 1 ^, modulate its specific heatil, affect 
the plasmon excitations 1 ^ and exhibit rich interplay with 
electron- electron interactions 1 ^. 

For technological applications, it is also desirable to 
understand and control its magnetic properties besides 
tunable electronic and optical properties. In order to 
investigate the intrinsic magnetism in graphene, Peres 
et al*& analyzed the exchange instabilities induced by 
the long- range Coulomb interaction between the Dirac 
fermions within the Hartree- Fock (HF) approximation. 
Their results showed that in undoped graphene, the 
itinerant electrons exhibit a ferromagnetic transition 
when the coupling constant is sufficiently large. However 
according to recent experiments 21 , no evidence for 
intrinsic magnetism was found at any temperature down 
to 2 K. This might be a consequence of weak effective 
interaction^ 2 - between the Dirac fermions and is a subject 
of on- going research. 

Apart from ferromagnetic exchange instabilities 
due to long- range interactions, there are numerous 
studies on the induced localized magnetism due to 
on- site Hubbard- type interactions in finite sheets 
(nanoribbons2£) and with dopants (vacancies and exter- 
nal adatoms^ 4 -) in graphene. Uchoa et al*^ examined 
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the possibility of local moment formation on adatom 
(impurity ion) in graphene. Their analysis was based 
on Anderson's single impurity model^ in which the 
constant density of states of the impurity hybridizes 
with the sea of conduction Dirac electrons in graphene. 
They concluded that under certain conditions (physical 
parameters of the model) it is feasible to form a localized 
magnetic moment on the impurity. They also showed 
that the local magnetic moment can be tuned by apply- 
ing a potential through an electric field via back gate, 
emphasizing the importance of electrically- controlled 
magnetic properties of an impurity in graphene. 

The purpose of this work is to explore the effect of 
an uniaxial strain not only on the itinerant magnetism 
i.e., ferromagnetic exchange instability due to long- 
range Coulomb interaction but also on the localized 
magnetism i.e., formation of localized magnetic moment 
due to an impurity in graphene. There are studies which 
aimed towards understanding magnetism in strained 
graphene but they are based on finite size nanoribbons 27 , 
topological line defects^, transition- metal atoms ad- 
sorbed on graphene based on density functional theory 
calculations^ or modulation of magnetic (RKKY-type) 
interactions between localized moments in graphene^. 

The rest of the paper is structured as follows. In 
section |n] we present the tight- binding Hamiltonian 
in the presence of an uniaxial strain and derive the 
anisotropic electronic dispersion of graphene. In sec- 
tion IIIII we examine the possibility of ferromagnetic 
exchange instabilities due to long- range Coulomb 
interaction between anisotropic Dirac fermions under 
applied strain within the HF approximation. We study 
the change in the strength of effective critical interaction 
or coupling constant which causes the magnetic transi- 
tion. In section IIVI we formulate the conditions for the 
existence of localized magnetic moments on an adatom 
in a deformed graphene. We present phase diagrams 
for different values of model parameters. In the last 
section [V] we summarize and conclude our findings. 



II. ANISOTROPIC ELECTRONIC 
BANDSTRUCTURE 

We begin by deriving the expression for the electronic 
dispersion of deformed graphene. In the undeformcd 
case, P. R. Wallace^ was the first to study its band struc- 
ture, as early as 1947, in terms of tight- binding model. 
In graphene the carbon atoms are arranged in a two di- 
mensional (2D) honeycomb lattice, see Fig. [2(a), with 
two atoms (A and B) in the unit cell. The Hamiltonian 
as obtained from tight- binding model reads 

H o = ^2 *(7i) a R<7 b R-+7i,<7 + H.c (1) 

where t("/i) = ti are the nearest neighbor hopping pa- 
rameters as shown in Fig. [Ha). In the isotropic case 
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FIG. 1: (Color online) Left : Undeformed graphene honey- 
comb lattice structure showing two sub- lattices A and B with 
ti, t2 and t3 as the nearest neighbor hopping parameters. 
Right : Uniaxially strained along the armchair, y— direction. 

they are all equalivalent ; to = U f=s 2.7eV. Here aj^ 
(bR+ 7ijCT ) is creation (annihilation) operator of an elec- 
tron at sub- lattice A (B) with spin a =t,| at posi- 
tion R (R + 7i) and defined as a^ CT — J^k e ~' k Ra Lr 
(b R+7lff = jf Ek elk (R+7l)b k<r)- The bond vectors from 
atom A (or B) to its nearest neighbor B (or A) are 

denoted by 7, with 71 = a(^,^), 72 = a(0, 1) and 
73 = a( = ^-, ^j.). The bond length or the distance be- 
tween two carbon atoms is a(«0.142 nm) and a = \J~Sa 
as shown in Fig. [IJa) . We consider uniform bond defor- 
mations and neglect the bond bending effects. 

In the momentum space the Hamiltonian in Eq. ([1]) 
can be written as 

H Q = -ta ^( k )4«A- + H c (2) 

where </>(k) = e lk ' 72 [l + 2rje~ lk y a y cos(k x a x )], 77 = |^ = 

ffj a a; = § an d ci y = ^a. We obtain the low- en- 
ergy effective Hamiltonian by expanding the wave vec- 
tor, k, near one of the six Dirac points. In fact there 
are only two, out of six Dirac points, which are inequiv- 
alent. Therefore on expanding around either of the two 
inequivalent Dirac points, (±^t=,Q) we get 

H = V X k X T X + VykyTy (3) 

where f x , f y are the usual (2 x 2) Pauli matrices. Here 
v x = t 2 a x ^A7f — 1, v y — tidy are the velocities along 
the x— and y— directions respectively. And on diag- 
onalizing the Hamiltonian in Eq. ([3]), we get the re- 
quired anisotropic electronic dispersion (elliptical double 
cones )i£ 

e(k) = ±yjvlk*+vffi (4) 

where ± denotes the electron- and hole- like bands. We 
define momentum cut- offs along x— and y— directions 
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in such a way that the number of states in the Brillouin 
zone are conserved while going from isotropic to the 
anisotropic case i.e., the area of the unit cell remains con- 
stant. Thus A£ = ^ = x f^-k c and M=^ = .f^k c 

c Vx y «i L C Vy V V V 

where D a is the anisotropic energy cut- off and k c being 
the isotropic momentum cut- off. And the area, A, of 

the unit cell is defined so that ^ 2n } = -^f- = k^M = k 2 ,. 

In the following sections, we inspect the conditions 
required to observe an exchange instability due to 
long- range Coulomb interaction (section IIIip and the 
possibility of forming a localized magnetic moment on 
an impurity (section IIV|) in deformed graphene. 



III. FERROMAGNETIC EXCHANGE 
INSTABILITY 

In 1929, Bloch proposed the occurence of ferromag- 
netic transition driven by electronic density in a three- 
dimensional (3D) interacting electron gas using Hartree- 
Fock mean- field theory 3 ^. The result was that at high 
density the electron system would have paramagnetic 
ground state in order to optimize the cost in kinetic en- 
ergy due to large number of electrons while at low den- 
sity the system should spontaneously spin- polarize into 
a ferromagnetic ground state so as to compromise for 
the exchange energy arising from the Pauli principle and 
Coulomb interaction. 

For an electron gas in a positive background, it is easy 
to calculate 3 - 3 , the total HF energy per particle (which 
is a sum of the non- interacting kinetic energy and the 
Fock exchange energy due to unscreened Coulomb inter- 
action at zero temperature) and see that it leads to first- 
order ferromagnetic transition, so- called Bloch transi- 
tion. This happens when the dimensionless parameter, 
r s which is defined as the ratio of potential to the kinetic 
energy and which is inversely proportional to the elec- 
tronic density, becomes greater than a critical value. In 
such a situation, the ferromagnetic (fully spin- polarized) 
state is lower in energy than the paramagnetic (unpolar- 
ized) state. 

In 2005, Peres et al*° investigated the possibility of 
such a transition in pure (half- filled undoped) graphene 
within the Hartrcc- Fock approximation. They concluded 
that for coupling constant, a, greater than certain critical 
value, a c ~ 5.3, the ground state of the system becomes 
maximally spin- polarized. Since the behavior of elec- 
trons in graphene is quite different than the usual two- 
dimensional electron gas, one might argue that it is in- 
adequate to calculate the ferromagnetic exchange insta- 
bility based on the Hartree- Fock approximation and one 
has to include the correlation effects i.e., going beyond 
the Hartree- Fock approximation. Such an estimate of 
the correlation energy was made by Dharma-wardana 34 , 
The author found out that the inclusion of correlations 
suppressed the exchange- driven spin- polarized phase. 



But undoped graphene has vanishing density of states 
near the Fermi energy (Dirac point) and the long- range 
Coulomb interaction remains unscreened. Moreover un- 
like usual electron gas system, the coupling constant in 
graphene is independent of its electronic density. So it is 
an elusive task to quantitatively establish the correlation 
effects in graphene. Following Ref. [2(], we examine the 
effect of uniaxial strain on the ferromagnetic exchange 
instabilities in graphene within Hartree- Fock approxi- 
mation. Even though it must be kept in mind that at 
Hartree-Fock level the tendency toward itinerant mag- 
netism is possibly too strong, as is the case in the conven- 
tional electron gas^ 3 ., this level of approximation allows 
for complete analysis of the effects due to strain (Dirac 
fermion anisotropy). 

We consider the case of charge neutral graphene, i.e. 
zero chemical potential. It is expected that, as in un- 
strained graphene^, the ferromagnetic instability is the 
strongest in this case. The ground state energy per par- 
ticle, per valley, of a spin polarized system is expressed in 
terms of the kinetic, exchange and correlation energies. 
The kinetic energy, from Eq. ([J}, is given by 

K = g s hY,\Jv 2 x k 2 x + vlkl 

k 

= 9s hA J I ^2 fc 2 +u 2 fc 2_|*_ (5) 

where g s —2 is spin degeneracy factor. It is straight- 
forward to derive the kinetic energy of the unpolarized 
(paramagnetic) state, K pi as a function of isotropic mo- 
mentum cut- off, k c ; 

K p = -^r^-V v x v y k c ( 6 ) 

and also the kinetic energy of the polarized (ferromag- 
netic) state, Kf, for either of the spin as a function of k c 
and the Fermi wavevector, fcj? ; 

HA 

Kf = -—^V^{kl-k F ) (7) 

Here the Fermi wavevector k F is related to the magne- 
tization M via: M = = (|^) 2 with n a being the 
occupation number for spin a . The kinetic energy differ- 
ence, AK , between the ferromagnetic and paramagnetic 
phases is given as 

AK = g s K f - K p 

fJsflA n , . 

= —„ y/VxVykp (8) 

07T 

where the spin degeneracy factor takes into account the 
individual contributions of either of the spins for the spin- 
polarized case. 

We now evaluate the interaction energy which consists 
of both exchange and the correlation terms. But as men- 
tioned earlier, we shall calculate only the exchange energy 
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term and neglect the correlation energy. In graphene, the 
electrons interact via long- range Coulomb interaction 

9 2 

which is defined as V p — where k is the appropriate 
dielectric constant. The exchange energy is determined 
by electron- hole excitations and within the Hartree- Fock 
approximation is given by 

k,p,tr 

(9) 

where G"(k, lu) is the non- interacting Green's function 
for spin a which is defined as 



G°(k,u;) = 



WTO + V x k X T x + Vytyfy 



- (v%k% + v*k$) 



(10) 



with fo being the (2 x 2) identity matrix and f x , r y are the 
Pauli matrices. The momentum summations in Eq. (|9]) 
are reduced to integrals with limits of the integrations 
as k c or &f depending on whether the exchange energy 
is calculated for paramagnetic or ferromagnetic phase. 
After making use of scaling transformations we evaluate 
the exchange energy difference, AE ex , between the fer- 
romagnetic and paramagnetic phases as 



AE r , 



hA^v x v y k F a x 



16tt 3 



h 



2/ 2 
v/M 



(11) 



where a x = - £ — is the strength of the coupling constant 
which is also given as the ratio of potential to the kinetic 
energy along x— direction with the directional depen- 
dence coming because of anisotropy, and 



h = 

and 

h 



J 



27T r 2-K 



rt(l + cos(6» - <f>))drdtd6d<p 



rtcos(6 — <fi)drdtd6d(f) 



io Jo Jo jo g(r,t,0,(f)) 
arc the integrals which depend on the functions 



(12) 



(13) 



f(r,t,0, </>)-- 
and 

g{r,t,9 i( f>) 



\ I —(rcosO — tcos(f)) ) + (rsin6> — tsm<j)) 2 
V \ v * J 



— (ry/Mcosd — teos0) ) + (ry/Msmd - tsin(/)) 2 
K v x J 



In what follows we shall define the anisotropy param- 
eter due to applied uniaxial strain. We define it to be 
proportional to the ratio of the two velocities, v y and 
v x i.e., let — = 1 — S. We are interested in the range 
< 5 < 1, i.e. 1 > v y /v x > 0, which according to 



Ferromagnetic 



or 



Paramagnetic 



1 - (v /v ) 



FIG. 2: (Color online) Boundary separating the paramagnetic 
and ferromagnetic regions indicating the onset of an exchange 
instability in deformed graphene. The value of the critical 
coupling constant, a%, decreases as a function of applied uni- 
axial strain. 



Fig. [ljb) suggests that the tension is along the armchair, 
y— direction with v y < v x . The case S — corresponds 
to the isotropic and S = 1 is the limit of decoupled chains 
of carbon atoms. In the similar manner we can describe 
strain along zig- zag, x— direction, by using the same 
parameter range and accordingly relabeling the axes. 

The total energy gain between the ferromagnetic and 
paramagnetic phase is given by the sum of the kinetic, 
Eq. (JSJ), and exchange energy differences, Eq. (fTTj); 



AE = AK + AE ex 
= hAy/v x v y k 3 F 



1 

37 



16tt 3 



2/ 2 



(14) 



where the integrals are functions of anisotropy and mag- 
netization. 

We find numerically that the exchange energy exceeds 
the kinetic energy for large magnetization M = 1 (i.e. 
kp = k c ), which leads to the interaction-driven ferro- 
magnetic transition. The case of partial spin polarization 
M < 1 is always energetically less favorable. The large 
magnetization is characteristic feature of the Hartree- 
Fock approximatio n 20 ' 33 . Such a scenario occurs when 
the coupling constant strength, a x , is greater than cer- 
tain critical value, ol% : 



Gt-f) Girt* 



16tt 2 



(15) 



Thus when the condition in Eq. ()T5|) is satisfied for 
a given anisotropy, the system exhibits a first-order 
transition (Bloch transition) and the spin degeneracy 
is lifted, i.e. for a x < a x the system is unpolarized 
(paramagnetic) and when a x > a c x it is fully spin- 
polarized (ferromagnetic) . 

We solve for the critical values of coupling constant 
for M = 1 and plot it as a function of anisotropy, 
8 = 1 — as shown in Fig. [5] The phase boundary, 



where the Bloch instability occurs, separates the unpo- 
larized from the fully spin- polarized region. Our main 
result is that with an applied uniaxial strain (increasing 
anisotropy) the value of the critical coupling constant, 
a%, decreases. In the limit of extreme strain, 6 = 1, the 
two- dimensional graphene decouples into chains of car- 
bon atoms. It is found that the value of critical coupling 
constant is finite and is depicted as a circle in Fig. [3] 
The reason a% decreases with increasing anisotropy is 
due to the fact that electronic instabilities are generally 
enhanced in lower dimensions. However we should still 
keep in mind that even for large anisotropy the critical 
interaction strength remains relatively large a% < 2 
and thus correlation energy effects could be substantial. 
Nevertheless our result demonstrates that itinerant 
ferromagnetism is strongly favored in strained graphene 
- an exciting possibility not yet tested experimentally. 
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FIG. 3: (Color online) Schematic of (a) an uniaxially strained 
graphene with an adatom (in green) sitting on top of one of 
the two sub- lattices (red or blue) and (b) the corresponding 
anisotropic dispersion of the Dirac fermions near one of the 
Dirac points along with flat dispersionless level of the localized 
adatom which can be moved in the conduction (as shown in 
the figure) or in the valence band. 



IV. LOCALIZED MAGNETIC STATES ON 
ADATOM 

In this section we investigate the conditions required 
to form localized magnetic moments on an adatom in 
strained graphene. Our work closely follows the analysis 
of Uchoa et a/. 2 — in terms of the description of the model 
using the single impurity Anderson Hamiltonian. Al- 
though their investigation was extended for the situation 
with magnetic adatom 35 and for an impurity on bilayer 
graphene 3 ^, the case of deformed monolayer graphene 
with a single non- magnetic adatom remains unexplored. 

We consider an uniaxially strained graphene sheet with 
a single adatom (impurity ion shown in green) sitting on 
top of one of the two incquivalcnt atoms (A or B shown 
in red or blue respectively) as shown in Fig. . The 
anisotropic dispersion (elliptical cone) for one of the val- 
ley degrees of freedom along with flat dispersionless level 
of the localized impurity in the conduction band is pre- 
sented in Fig. [3jb) . The model Hamiltonian for such a 
system in momentum space is expressed as 



H — Hq + Hf + H\ 



(16) 



where the first term , Hq, is the anisotropic dispersion of 
the Dirac electrons in graphene as given in Eq. (|2|) . The 
second term describes the featureless localized orbital of 
an impurity ion (constant density of the states at energy 
i?o) along with the Coulomb repulsion energy which takes 
into account the double occupancy of the localized energy 
level, 



level, nl, depends on its density of states and will be 
defined in the following text. The localized level is pop- 
ulated by the Dirac electrons due to hybridization which 
is given by the last term in Eq. JTF 



Hv = V^2(fta pa 

pa 



"pa 



fa) 



(18) 



where V is the strength of hybridization and we assume 
that the impurity sits on top of atom A. 

We use the Green's function method to solve the 
Hamiltonian in Eq. (1161) . Since we are interested in the 
condition of formation of localized magnetic moment on 
adatom, we consider the full interacting Green's function 
of the electrons in the localized level which is defined as 
G fa {uj) = [cj-e CT -£ / (w)+iO + ]~ 1 where e a = E + Unt a 
is the energy of the localized electrons with spin a in the 
presence of interaction energy U and the self-energy has 
the form. 



lo 2 - [(v x p x ) 2 + (v y p y ) 2 



D 2 



udoe 



D 2 



t7r|w|e(D a - M) 

(19) 



(E + UnUfU 



(17) 



where /J creates (annihilates) an electron with spin 
a at the impurity site and we have used the mean- field 
decoupling of the local interaction 25 , U. The spin- depen- 
dent occupation number of the electrons at the impurity 



where D a is the anisotropic energy cut- off as defined in 
section ITlland O being the Heaviside step function. With 
the knowledge of the full Green's function we can obtain 
the density of states of the impurity electrons in the lo- 
calized state, Pfai^) = — — Gf a (u)). 

The condition to form a magnetic moment on the im- 
purity is determined by calculating the spin- dependent 
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FIG. 4: (Color online) Phase diagram showing the magnetic regions signaling the onset of formation of localized magnetic 
states on the adatom (impurity) in graphene for various values of applied uniaxial strain. The x— and the y— axis of the 
phase diagram are scaled model parameters. Each phase plot is for a fixed value of localized impurity energy level (-Eo) and 
the strength of bare hybridization (V) in units of eV. The plots in the upper (lower) panels are for positive (negative) values 
of So. 



occupation number which is given by 

Pf a {L0)d,L0 



nl = 



A /-k 

7T 



\x\dx 



1 - #log(T^) 



A 2 x 2 



(20) 



where Ep is the Fermi energy which is tunable param- 
eter by applying gate voltage, A = is the effective 

hybridization and £ a — Eo+ ^ n -" j s ^he energy of the 
localized electrons which is scaled with respect to the 
anisotropic cut- off. It is evident from Eq. (|20[) that the 
determination of spin- dependent occupation number re- 
quires the self- consistent calculation of the density of 
states, pf a (uS), which incorporates the broadening of the 
impurity level due to hybridization with the sea of Dirac 
electrons in graphene. A localized moment is formed at 
the impurity whenever there is an imbalance in the spin- 



dependent occupation number at the level. 

We define anisotropy parameter due to strain through 
the anisotropic energy cut- off, D a which depends on 
velocities; v x and v y . When we apply tension along 
the armchair (y— ) direction as shown in Fig. [IJb) the 
bond vectors yields |7i| = I73I = (1 + | — ^-)a and 
I72I = (1 + e)a where we assume that the change in 
bond length is linear in strain^ (e) and v — 0.165 is the 
Poisson's ratio of graphite^. On taking the functional 

form of hopping parameters^ as ti — t^e Aa ( a / , we 
get the velocities as a function of applied strain which 
are given by v x = vfO- + Xaeu) and v y — vp(l — Xas) 
where vf — 2^ is the isotropic Fermi veolcity. The 
anisotropic energy cut- off as a function of strain is given 
by, D a = y/v x V y k c w D(l — ce) with D = vpk c m 7 eV be- 
ing the isotropic graphene bandwidth and c = Aa(l— v) is 
a constant with Xa ps 3.37 as given in Ref.fl2l Notice that 
the anisotropy parameter is defined through the product 
of two velocities, v x and v y , unlike in the previous section 
where it was a function of ratio of the two velocities. 
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We regard x = -jj- and y = ( F ^ 0> as scaled model 
parameters and solve Eq. (50] for different values of im- 
purity level (-Eq), strength of bare hybridization (V) and 
local interaction (U) given in units of eV and for applied 
strain (e). We obtain a phase diagram as shown in Fig.0] 
and observe clear boundary separating the magnetic and 
non- magnetic regions. Within our model the impurity 
level can be shifted either in the conduction band i.e., 
Eq > as shown in Fig. EJb) or in the valence band i.e., 
Eq < and the corresponding phase plots are shown in 
the upper or lower panels of Fig. [4] respectively. 

For the sake of completeness, we briefly summarize the 
main results of undeformed graphene (e = 0) as obtained 
in Ref. [25]. As seen in Fig. [4j the phase diagram is not 
symmetric for the cases where the localized energy level 
is either above or below the Dirac point and also around 
y = 0.5 which shows the particle- hole symmetry break- 
ing due to existence of the impurity level. Moreover, 
when Eq > the magnetic boundary crosses the line 
y = and the impurity magnetizes even when the en- 
ergy level is above the Fermi energy. And similarly for 
the case E < 0, the crossing occurs along y = 1 which 
suggests that even when the energy of the doubly oc- 
cupied level is less than the Fermi energy the impurity 
gets magnetized. This happens due to the fact that the 
hybridization leads to large broadening of localized level 
density of states that crosses the Fermi energy. These 
features are very different as compared to the case of an 
impurity in a metal 25 . In the following, we systemati- 
cally study the effect of uniaxial strain on the the scaling 
of the magnetic phase boundary. 

We first consider the case when the localized level is 
in the conduction band (Eq > 0) and show the results 
in Fig. U{a) - ( c )- For a fixed strength of bare hybridiza- 
tion (V — 0.2) shown in Fig. |U[a), which is considered 
to be weak as compared to the graphene bandwidth, the 
magnetic region starts depleting as the applied strain in- 
creases. And this effect is more pronounced with stronger 
hybridization (V — 1.0) as shown in Fig. B|b) and (c). 
The reason being with increasing strain the anisotropic 
cut- off energy (D a ) decreases and for given value of V, 
the effective hybridization (A) increases. Thus the den- 
sity of states at any given energy increases since it is pro- 
portional to the effective hybridization and this results in 
the electrons at the impurity level becoming delocalized. 
Moreover the derealization increases upon enhancing the 
hybridization strength, V. 

It is also observed that for undeformed case, the size of 
the magnetic region gets enlarged as the impurity level 
(Eq) approaches the zero- energy Dirac point for a given 
fixed hybridization. The formation of local moment is 
favored due to suppresion of density of states at low en- 
ergies. If the value of bare hybridization is large then 
with strong deformations the magnetic region gets dras- 
tically reduced, as seen in Fig. 0Jb) and (c), due to aug- 
mented density of states. In such a situation there is an 
emergence of a critical value of the scaled parameter (x) 
above which the impurity gets magnetized. This happens 



because the model breaks down for small x and substan- 
tial anisotropy, since in this case the Hubbard-shiftcd 
impurity level Eq + U is actually outside the Dirac cone 
bandwidth. 

The depletion in the magnetic region with applied 
strain is also found for the case when the localized level is 
in the valence band (Eq < 0) and the results are shown in 
Fig. Bid) - (f). Apart from similar observations as in the 
case of Eq > , it is also seen that for large hybridization 
and small anisotropy, there's an up turn close to point 
y = 1 and x m i.e., the magnetic phase boundary ap- 
proaches that point from below. This suggests that in 
such a scenario the system behaves more like an impu- 
rity in a metal, i.e. the usual single impurity Anderson 
model is regained. 



V. CONCLUSION 

In this paper we examined the effect of applied uniax- 
ial strain on the ferromagnetic exchange instability due 
to long- range Coulomb interaction between the Dirac 
fcrmions and also on the formation of localized magnetic 
moment on an adatom (impurity) in graphene. Our work 
was motivated by a broad aim to understand the behavior 
of itinerant as well as localized magnetism in deformed 
graphene. 

We derived the electronic dispersion of the anisotropic 
Dirac fermions and considered that the electrons in 
graphene are interacting via long- range Coulomb inter- 
action. We calculated the total energy difference between 
the ferromagnetic (polarized) and paramagnetic (unpo- 
larized) phases within the Hartree-Fock approximation 
and determined the necessary condition for the occurence 
of the exchange instability or the paramagnetic to ferro- 
magnetic transition as a function of anisotropy. In the 
undeformed case, Peres et al. conjectured that large criti- 
cal coupling constant, i.e. strength of interaction can lead 
to first order magnetic (Bloch) transition. We found that 
the value of critical coupling constant, which caused ex- 
change instability, decreased with increasing applied uni- 
axial strain. For large deformations our work indicates 
that graphene can have a ferromagnetic ground state for 
relatively weak effective interaction (a < 2) which is com- 
parable to the case of suspended graphene^,. The effect 
of uniaxial strain on the electronic properties has already 
been studied experimentally but the magnetic properties 
of deformed monolayer graphene have not been explored. 
It would be very interesting to confirm our result exper- 
imentally. 

We also formulated the conditions for the existence of 
localized magnetic moments on an adatom in a deformed 
graphene. Using the single impurity Anderson model and 
within mean- field self- consistent approach, we system- 
atically studied the polarization of the electrons in the 
localized level as a function of model parameters like hy- 
bridization strength (V), on- site Coulomb repulsion (U), 
localized energy level (Eq) and the applied strain (e). We 
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obtained magnetic/non- magnetic phase boundary for dif- 
ferent values of scaled model parameters and found that 
for moderate hybridization the magnetic region shrunk 
drastically with increasing strain as compared to weak 
hybridization. This can be understood by the fact that 
the anisotropic cut- off energy kept on decreasing with 
increase in deformation, which resulted in an increase in 
effective hybridization (A). This in turn increased the 
density of states at the localized level and the electrons 
at the impurity level became more de-localized. The fea- 
tures predicted in this work can find possible applications 



in the field of carbon-based spintronics. 
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